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Better  Pasting  via  Quasi-Interpolation 

Blair  Conrad  and  Stephen  Mann 


Abstract.  Surface  pasting  is  a hierarchical  modeling  technique  that 
places  feature  surface  atop  base  surface  with  approximate  continuity.  In 
this  paper,  we  describe  the  use  of  quasi-interpolation  to  set  the  bound- 
ary control  points  of  the  pasted  feature.  To  achieve  interactive  perfor- 
mance for  surface  pasting,  we  modify  the  coefficients  and  weights  of  quasi- 
interpolation and  devise  an  efficient  sampling  scheme.  The  new  sur- 
face pasting  technique  is  actually  faster  than  the  original  surface  pasting 
method  when  using  the  same  number  of  feature  control  points,  with  a 
much  lower  discontinuity  between  the  feature  and  base  surface. 


§1.  Introduction 

Spline  curves  and  surfaces  are  used  in  many  areas  of  computer  graphics  and 
computer  aided  geometric  design.  In  particular,  tensor  product  B-spline  sur- 
faces are  commonly  used  in  modeling  and  computer  animation  because  they 
have  many  attractive  properties,  such  as  a compact  representation  and  ad- 
justable levels  of  internal  continuity  [7]. 

Frequently,  the  user  of  a piece  of  modeling  or  animation  software  will 
want  to  add  a region  of  local  detail  to  a B-spline  tensor  product  surface, 
but  the  knot  structure  will  be  too  coarse  to  allow  the  fine-grained  control 
that  the  user  desires.  Traditional  methods  of  increasing  the  complexity  of 
the  surface  include  inserting  knots  using  either  Boehm’s  algorithm  [3]  or  the 
Oslo  algorithm  [6].  The  insertion  of  a knot  into  either  of  a surface’s  knot 
vectors  causes  an  entire  row  or  column  of  subpatches  to  be  split  — rather 
than  increasing  the  number  of  subpatches  locally,  extra  subpatches  are  created 
across  the  width  or  breadth  of  the  surface. 

Forsey  and  Bartels  [8]  developed  hierarchical  B-splines  to  allow  the  user 
to  add  local  detail  to  a tensor  product  B-spline  surface  while  maintaining  B- 
spline  continuity.  Hierarchical  B-splines  suffer  from  several  drawbacks:  the 
regions  of  added  detail  must  remain  parametrically  aligned  with  the  base, 
and  it  is  impractical  to  slide  the  features  or  to  maintain  a library  of  overlays 
that  may  be  added  to  a base  surface.  Wavelets  [10]  can  be  used  in  a similar 

Curve  and  Surface  Design:  Saint-Malo  1999  27 

Pierre-Jean  Laurent,  Paul  Sablonnifere,  and  Larry  L.  Schumaker  (eds.),  pp.  27-36. 

Copyright  © 2000  by  Vanderbilt  University  Press,  Nashville,  TN. 

ISBN  0-8265-1356-5. 

All  rights  of  reproduction  in  any  form  reserved. 


28 


B.  Conrad  and  S.  Mann 


hierarchical  modeling  technique,  but  they  also  require  the  detail  regions  to  be 
parametrically  aligned  with  the  base. 

Displacement  mapping  is  another  technique  for  adding  local  detail  to  a 
surface  while  adding  as  few  extra  control  points  as  possible.  Every  point  on 
the  feature  is  defined  by  a displacement  vector  relative  to  a certain  point  in 
a reference  plane.  To  evaluate  the  surface  at  the  detail,  you  add  the  displace- 
ment vector  of  the  feature  to  the  corresponding  point  on  the  base  surface. 
The  greatest  disadvantage  to  using  displacement  mapped  surfaces  is  the  com- 
putational cost.  Multiple  surface  evaluations  must  be  used  to  determine  each 
point  to  be  rendered  on  the  composite  surface. 

Surface  pasting,  suggested  by  Forsey  and  Bartels  [2],  is  a generalization  of 
hierarchical  B-splines  that  was  intended  to  combine  the  flexibility  of  displace- 
ment mapped  surfaces  with  the  speed  of  evaluation  enjoyed  by  hierarchical 
B-splines.  This  technique  has  the  flexibility  of  displacement  mapping,  but  it 
is  much  cheaper  since  only  the  control  points  of  the  feature  must  be  mapped, 
rather  than  the  larger  number  of  surface  points  to  be  rendered.  This  combina- 
tion of  flexibility  and  speed  has  drawn  the  attention  of  the  modeling  industry 
to  surface  pasting  — recent  versions  of  Houdini,  a commercial  animation  tool 
produced  by  Side  Effects  Software,  have  included  support  for  surface  pasting. 

However,  surface  pasting  is  only  an  approximation,  and  as  such  it  does 
not  have  the  same  continuity  properties  as  displacement  mapping  or  hierar- 
chical B-splines.  In  general,  there  is  no  guaranteed  continuity  between  the 
feature  and  the  base  surfaces.  By  inserting  knots  into  the  feature  surface, 
the  error  between  the  feature  boundary  and  the  base  surface  can  be  reduced 
to  any  desired  tolerance,  but  many  knot  insertions  may  be  required  to  get 
the  desired  approximation,  and  the  resulting  additional  control  points  in  the 
feature  dramatically  increase  the  cost  of  the  pasting  operation. 

In  this  paper,  we  suggest  altering  the  surface  pasting  technique  to  improve 
the  approximate  continuity  between  the  feature  and  base  surfaces  by  using 
quasi-interpolation  to  set  one  or  more  of  the  outer  rings  of  feature  control 
points. 


§2.  Surface  Pasting 

Since  the  work  in  this  paper  is  an  extension  of  surface  pasting,  we  will  give  a 
quick  review  of  how  pasting  works.  For  details  on  standard  surface  pasting, 
see  any  of  the  earlier  works  on  the  subject  [1,2]. 

The  pasting  process  is  illustrated  in  Figure  1.  Surface  pasting  combines  a 
base  surface  and  a feature  surface,  each  of  which  is  in  tensor-product  B-spline 
form.  The  feature’s  control  points  are  adjusted  so  that  the  boundary  of  the 
pasted  feature  lies  on  or  near  the  base  surface,  and  the  shape  of  the  pasted 
feature  reflects  the  original  shape  of  the  feature  imposed  as  a displacement  of 
the  base  surface. 

To  map  the  feature’s  control  points,  the  feature’s  domain  is  embedded  in 
the  feature’s  range  (upper  left  of  Figure  1).  Next,  a local  coordinate  frame 
= { ui,j  i vi,j , wi,j  i Oij  } is  constructed  for  each  feature  control  point 
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Fig.  1.  Surface  Pasting. 

where  the  origin  Oij  of  each  frame  is  the  Greville  point  corresponding  to 
PiJ,  two  of  the  frame’s  basis  vectors  are  the  parametric  domain  directions 
and  the  third  basis  vector  is  the  direction  perpendicular  to  the  domain.  Each 
control  point  P2tJ  is  then  expressed  relative  to  its  local  coordinate  frame 
as  Pii:j  = auij  + Pvij  + yiVij  + Oi}j. 

Next,  we  associate  the  feature’s  domain  with  a region  in  the  base’s  domain 
(right  half  of  Figure  1).  This  gives  the  location  on  the  base  surface  where  we 
will  place  the  feature.  We  now  map  each  coordinate  frame  Tij  onto  the  base 
surface,  giving  a new  coordinate  frame  P-j  = whose 

origin  O'ij  is  the  evaluation  of  the  base  surface  at  Oi,j,  and  two  of  its  basis 
vectors  lie  in  the  tangent  plane  of  the  base  surface  at  that  point,  the  third 
being  perpendicular  to  the  tangent  plane.  We  then  use  the  coordinates  of 
each  feature  control  point  Pij  relative  to  as  coordinates  of  the  elements 
of  the  frame  T[y  This  gives  us  the  location  of  the  pasted  feature  control 
point,  PP  = au'itj  + Pv[j  + 7 w'id  + O'ij. 

§3.  Feature  Boundaries 

With  standard  surface  pasting,  the  boundary  control  points  of  the  pasted  fea- 
ture surface  (Figure  2)  are  conventionally  placed  on  the  base  surface,  resulting 
in  an  almost  C°  join.  Likewise,  the  second  layer  of  feature  control  points  are 
also  conventionally  located  on  the  base  surface,  giving  an  almost  C1  join.  By 
inserting  knots  into  the  feature  surface,  the  discontinuity  between  the  feature 
and  the  base  can  be  made  as  small  as  desired. 

One  disadvantage  of  using  knot  insertion  to  decrease  the  discontinuity 
between  the  feature  surface  and  the  base  surface  is  that  it  greatly  increases  the 
number  of  control  points  in  the  feature  surface.  In  addition  to  increasing  the 
storage  requirements  of  the  feature  surface,  this  increases  the  computational 
cost  for  interactive  modeling,  as  each  feature  control  point  needs  to  be  mapped 
with  the  pasting  process  described  in  the  previous  section. 

The  goal  of  this  work  is  to  find  better  settings  of  the  boundary  layers  of  the 
feature  control  points  to  minimize  the  C°  and  C1  discontinuities.  Thus,  we  will 
find  better  settings  of  the  control  points  illustrated  in  Figure  2.  In  this  figure, 
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Fig.  2.  Approximate  C°  boundary  points  and  C1  boundary  points. 

the  circled  points  are  set  using  quasi- interpolation:  the  black  control  points 
will  be  set  to  interpolate  position  and  derivative  information  sampled  from 
the  base  surface;  the  gray  control  points  will  be  set  using  quasi-interpolation; 
and  the  remaining  control  points  will  be  set  using  the  standard  surface  pasting 
process. 


§4.  Quasi-Interpolation 

Quasi-interpolation  is  a spline  approximation  technique  developed  by  de  Boor 
and  Fix  [4],  Their  method,  given  a function  / defined  over  a region  of  7?. 
and  a partition  tt  of  TZ,  constructs  a degree  m spline,  F^f,  that  approximates 
/.  FV/  is  called  the  quasi-interpolant  of  /.  The  quasi-interpolant  is  a local 
approximation  in  that  its  value  at  u depends  only  on  the  values  of  / in  a small 
neighborhood  around  u,  it  reproduces  polynomials  of  degree  m or  less,  and  it 
provides  a high  order  approximation  to  /,  with  \F„f  - f \ being  <3(Mm+1). 

Lyche  and  Schumaker  [9]  describe  a family  of  quasi-interpolation  opera- 
tors Q which  include  the  one  of  de  Boor  and  Fix.  They  take 

M 

(1) 

1=0 

where  the  {B™}fL0  are  the  B-splines  and  {A;},^.0  are  linear  functionals.  The 
linear  functionals  can  be  chosen  so  that  Q is  applicable  to  a wide  class  of 
functions,  Q is  local,  and  Qf  approximates  smooth  functions  with  a high  order 
of  accuracy.  Then  the  operator  Q can  be  applied  to  a real-valued  function  / 
to  produce  a B-spline  curve  Qf  that  approximates  /. 

The  quasi-interpolant  Qf  approximates  a smooth  function  / with  a high 
order  of  accuracy  because  Q is  constructed  specifically  to  reproduce  polynomi- 
als. The  Lyche-Schumaker  quasi-interpolants  approximate  smooth  functions 
with  accuracy  of  up  to  0(hm+1)  for  a function  with  a sufficiently  high  degree 
of  continuity.  In  addition,  quasi-interpolants  can  be  constructed  using  linear 
functionals  that  can  be  quickly  evaluated,  resulting  in  a fast  approximation  op- 
erator. The  combination  of  high  accuracy  and  speed  make  quasi-interpolation 
a fitting  tool  to  improve  the  approximate  continuity  around  the  boundary  of 
pasted  features. 
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The  quasi-interpolant  Qf  is  in  B-spline  form.  Each  control  point  of  Qf 
is  a weighted  sum  of  linear  functionals  applied  to  /: 


Pi  — ^ ai,j^i,jf- 

l=o 

One  particular  choice  of  the  Lyche-Schumaker  quasi-interpolants  corresponds 
to 

^ i,j  f = 

where  the  r’s  are  locations  at  which  to  evaluate  /,  and  cxjj  is  the  blossom  [11] 
of 

Pi,j(u)  = (U  - Tii0)(«  - Ti} i)  ■ • ■ (U  - 

evaluated  at  iq, iq+m_i,  where  the  rq  are  the  knots  of  Qf.  This  results  in 
a quasi-interpolant  where  Qf  = f when  / is  a degree  m or  lower  polynomial 
and  if  Qf  ^ f,  then  the  error  has  the  best  possible  order. 


§5.  Modified  Quasi-interpolant 

The  Lyche-Schumaker  quasi-interpolant  uses  cheap  coefficients  and  expensive 
linear  functionals.  For  interactive  modeling  using  surface  pasting,  the  linear 
functionals  are  recalculated  frequently  (i.e.,  when  we  slide  the  feature  across 
the  base)  and  the  coefficients  are  recalculated  less  frequently  (only  when  we 
perform  knot  insertion  to  reduce  the  C°  discontinuity).  Thus,  to  match  the 
quasi-interpolant  to  our  application,  we  made  new,  cheaper  linear  functionals, 


^i,jf  — f{Ti,j)- 


This  choice  of  A ij  results  in  more  expensive  ctij,  which  are  now  the  blossom 
of 


Pi,j(u ) = n 

kjij 


U - Tj,k 
Ti,j  ~ Ti,k 


(2) 


evaluated  at  ?q, ..., where  the  are  the  knots  of  Qf.  Note  that  the 
Lyche-Schumaker  linear  functionals  require  computing  a divided  difference  of 
the  samples  while  ours  do  not,  but  the  denominator  term  of  (2)  makes  our  a 
coefficients  more  expensive  than  theirs. 

To  use  the  quasi-interpolant  on  all  four  boundaries  of  a feature,  the  ends 
of  the  connecting  boundary  curves  must  match.  This  is  most  easily  achieved 
by  modifying  Q to  reproduce  position  and  (to  use  the  quasi-interpolant  on  the 
cross-boundary  derivatives)  d derivatives  at  its  endpoints. 

Thus,  we  devised  a quasi-interpolant  that  we  call  Qd,  where  the  linear 
functionals  for  the  control  points  near  the  ends  of  the  curves  are  derivatives  of 
the  original  curve  at  the  endpoints.  This  quasi-interpolant  is  actually  a spe- 
cial case  of  the  Lyche-Schumaker  quasi-interpolant.  The  Qd  quasi-interpolant 
requires  different  A’s  and  p’s  at  the  endpoints,  which  in  turn  require  new  a’s. 
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For  0 < i < M,  and  0 < j < m,  we  define  Qd  as 


where 


M M m 

Qdf  - Ep*s.m  = E (Ea*.i\i/)5,m, 


1=0 

t=0  j= 0 

/Kj). 

if  d < i < M — d 

D3f(um-i 

),  if  i < d and  j < i 

D3 

if  i > M — d and  j < M — i 

[0,  otherwise 

( u~Titk 


Pij(u)  = < 


k^j 


Ti,j 


, if  d < i < M — d 


{u-um- i)j  . 

, if  i < a and  j < i 


(u-uMy 


if  i > M — d and  j < M — ; 


{0,  if  i < d and  j > i 

0,  if  i>  M — d and  j > M — i 

Pij(ui,...,Ui+m- 1),  otherwise 


(3) 


where  the  Tijs  are  as  defined  in  the  next  section.  Note  that  for  d < i < 
M — d,  the  pijs  are  merely  Lagrange  polynomials.  This  quasi-interpolant  Qd 
interpolates  the  position  and  derivatives  at  the  endpoints  and  has  optimal 
error  bounds  as  discussed  by  Conrad  [5]. 


Sampling  discipline 

The  above  equations  require  us  to  sample  the  function  / to  be  approximated 
at  some  values  Tij,  and  (3)  suggests  that  we  require  about  m + 1 samples 
of  the  base  surface  for  each  boundary  control  point  of  the  feature.  Since  the 
quasi-interpolant  error  bounds  hold  for  most  choices  of  r’s,  we  managed  to 
reduce  the  number  of  samples  of  the  base  surface  by  selecting  the  in  the 
following  way: 

• Place  the  first  d and  last  d control  points  in  groups  of  their  own;  divide 

the  remaining  control  points  into  groups  as  indicated  by  the  pseudo-code 
of  Figure  3;  each  group  k will  contain  control  points  Pffk , ...,  PNk+nk-u 
where  JVfc  = 11 1 and  is  the  number  of  control  points  in  group  k; 

• Let  7 i be  the  Greville  point  associated  with  control  point  Pt;  then  for 
each  group  k,  choose  the  interval  from  which  to  sample  as 

[0fcA]  = [(7Nt— i +7NJt)/2,(7wJt+1— i +7wt+,)/2]; 

• Set  Titj  = ak  + j{bk  ~ ak)/m,  for  i = Nk,  ...,Nk+i  - 1 and  j = 0, 

• Sample  / at  the  Tij\  note  that  you  can  share  the  first  and  last  samples 
between  groups, 
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I :=  M — 2 * d + 1 
numFullGroups  :=  [I /m\ 

L I — m * numFullGroups 
if  L = 0 

make  numFullGroups  groups  of  m 
else  if  numFullGroups  = 0 (mod  2) 

put  leftovers  in  middle  group;  all  others  have  m 
else  there  are  an  odd  number  of  full  groups 
if  L = 0 (mod  2) 

first  and  last  groups  contain  L/ 2 each;  all  others 

else 

move  one  control  vertex  from  middle  full  group  to  the  leftovers 

the  first  and  last  groups  each  contain  ( L + l)/2, 

the  middle  group  contains  m — 1,  and  the  others  contain  m 

Fig.  3.  Pseudo-code  for  grouping. 


Fig.  4.  Sampling  discipline  for  quasi-interpolated  surface  pasting. 

Pseudo-code  for  our  grouping  scheme  appears  in  Figure  3,  and  the  sampling 
is  illustrated  in  Figure  4.  The  white  and  gray  points  are  samples  of  the  base 
surface;  the  black  points  represent  the  control  points  for  the  quasi-interpolant 
(the  actual  quasi-interpolant  control  points  will  lie  much  closer  to  the  curve; 
we  have  moved  them  away  from  their  actual  position  for  illustrative  purposes); 
the  dotted  lines  indicate  the  boundaries  between  the  groups,  with  the  gray 
points  being  shared  by  two  groups;  and  the  triangles  indicate  the  position  of 
the  Greville  points. 

This  sampling  method  reduces  the  number  of  base  surface  evaluations  to 
about  one  sample  per  control  point. 

Quasi-interpolated  surface  pasting 

Ideally,  we  would  be  able  to  integrate  quasi-interpolation  into  surface  pasting 
by  running  four  independent  quasi-interpolations  to  set  the  boundary  control 
points,  then  four  more  independent  quasi-interpolations  to  set  the  second  layer 
of  control  points,  and  so  on.  However,  as  seen  in  Figure  2,  the  boundary  layers 
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of  control  points  overlap,  and  while  the  gray  control  points  in  this  figure  can 
be  set  independently,  the  black  ones  are  shared  by  two  boundaries.  To  build 
a tensor-product  patch,  we  need  the  two  quasi-interpolation  solutions  that 
compute  a group  of  black  points  to  give  these  points  the  same  settings. 

Our  solution  is  to  sample  the  base  surface  at  the  corners  of  the  embed- 
ded feature  domain  for  position  and  as  many  derivatives  as  we  are  trying  to 
approximate  along  the  boundary.  In  our  implementation,  we  were  only  trying 
to  achieve  approximate  C 1 continuity,  so  we  sampled  for  position  and  first 
derivatives  (including  mixed  partial  derivatives). 

Next,  we  used  the  Qd  operator  to  set  d + l outermost  rings  of  control 
points.  With  d = 1 (the  right  side  of  Figure  2),  this  means  that  we  make  eight 
applications  of  the  Q 1 operator,  four  for  the  outer  layer  of  control  points, 
and  four  for  the  next  layer  of  control  points.  In  each  application,  the  Q 1 
operator  interpolates  four  values  (those  indicated  by  the  black  points)  and 
quasi-interpolates  the  remaining  values. 

Note  that  while  we  quasi-interpolate  points  to  compute  the  first  layer  of 
control  points,  to  compute  the  second  layer  of  control  points,  we  first  quasi- 
interpolate  cross-boundary  derivative  vectors.  The  vector  coefficients  of  this 
second  quasi-interpolant  are  then  scaled  and  added  to  the  first  layer  control 
points  to  produce  the  second  layer  of  control  points. 

We  set  the  remaining  interior  control  points,  which  do  not  affect  the 
cross-boundary  derivatives,  using  the  standard  pasting  method. 

Results 

We  integrated  the  Q°  and  Q1  quasi-interpolants  into  surface  pasting.  In  our 
particular  implementation,  this  results  in  a cubic  quasi-interpolant  being  used 
to  approximate  a C2  curve.  Theoretically,  the  quasi-interpolant  described  in 
this  paper  should  have  0(h3)  order  of  accuracy,  which  matched  our  empirical 
tests.  This  is  an  improvement  over  standard  pasting,  which  empirically  had 
0(h2)  accuracy  on  the  same  data. 

We  achieved  a similar  improvement  for  the  cross-boundary  derivatives 
(improving  from  0(hl)  to  0(h2)),  and  the  Q1  method  has  an  additional  ad- 
vantage over  standard  pasting:  for  standard  pasting,  you  insert  knots  in  one 
parametric  direction  to  decrease  the  C°  discontinuity,  and  insert  knots  in  the 
other  parametric  direction  to  decrease  the  C1  discontinuity.  With  the  Ql 
method,  inserting  knots  in  one  parametric  direction  will  decrease  both  the  C° 
and  Cl  discontinuity  along  the  boundary. 

Visually,  the  quasi-interpolation  surface  pasting  method  gives  significant 
improvements  in  the  approximate  Cd  continuity  around  the  feature  bound- 
aries. An  example  appears  in  Figure  5.  In  all  three  images,  the  feature 
surface  has  the  same  number  of  control  points.  The  large  gap  that  appears  in 
the  standard  surface  pasting  example  has  almost  disappeared  with  the  Q°  and 
Q 1 pasted  surfaces,  but  a “corner”  has  appeared  in  the  Q°  surface.  However, 
Ql  pasted  surfaces  have  better  approximate  continuity  (both  C°  and  C1)  and 
are  cheaper  to  paste  than  standard  pasting,  assuming  we  use  feature  surfaces 
having  the  same  number  of  control  points  for  both  standard  and  Q1  pasting. 


Better  Pasting  via  Quasi-Interpolation 


35 


Fig.  5.  Standard  pasted  feature,  Q°  pasting,  and  Q1  pasting. 


This  computational  gain  for  Q1  pasting  is  a result  of  quasi-interpolating 
the  cross-boundary  vector  field  when  computing  the  second  layer  of  control 
points.  We  obtain  the  vectors  to  quasi-interpolate  when  we  evaluate  the  base 
surface  for  quasi-interpolating  the  boundary  curve.  Thus,  we  avoid  having  to 
evaluate  the  base  surface  when  computing  the  second  layer  of  control  points. 

The  biggest  computational  gain,  however,  comes  from  needing  fewer  con- 
trol points  for  Q 1 feature  surfaces  than  for  standard  pasting  surfaces  to  achieve 
the  same  tolerance  along  the  boundaries.  In  particular,  in  the  examples  we 
tested,  the  Q 1 method  required  one  third  the  number  of  boundary  points 
as  the  standard  surface  pasting  method  to  achieve  the  same  error  along  the 
boundary.  This  results  in  one-ninth  the  total  number  of  control  points  in  the 
feature,  resulting  in  roughly  a factor  of  ten  speedup  for  quasi-interpolated 
surface  pasting  over  standard  surface  pasting. 

Conclusion 

Surface  pasting  is  a flexible  modeling  technique  that  allows  the  creation  of 
multi-resolution  surfaces  by  hierarchically  composing  tensor  product  B-spline 
surfaces.  The  resulting  composite  surface  is  only  approximately  C°  and  ap- 
proximately C1.  In  this  paper,  we  developed  quasi-interpolation  operators  to 
reduce  the  C°  and  C1  discontinuities  between  the  feature  and  base  surfaces. 
By  reducing  the  C°  and  Cl  discontinuities,  we  can  use  fewer  control  points  in 
the  representation  of  the  feature  surfaces,  resulting  in  a significant  speedup  of 
the  pasting  technique. 

We  implemented  the  technique  described  in  this  paper  for  pasted  surfaces 
of  arbitrary  degree,  and  tested  the  method  on  cubic  and  quartic  surfaces. 
While  we  have  only  reduced  C°  and  Cl  errors,  the  extension  to  higher  levels 
of  approximate  continuity  is  straightforward,  but  calculating  the  coefficients 
that  weight  the  linear  functionals  will  be  more  complicated. 
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